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Abstract 

The  effect  of  cathode  particle  size  on  cell  potential  was  investigated  using  a  voltage  prediction  method  from  the  Coulomb  potential  created 
by  the  atoms  of  a  cathode  active  material.  LiNiC>2  was  selected  as  a  typical  cathode  of  an  ordered  rock  salt  structure.  From  the  calculation,  it 
was  found  that  the  voltage  changes  when  the  size  becomes  less  than  about  10  nm.  The  voltage  increased  with  decreasing  in-plane  size,  and 
decreased  with  decreasing  the  number  of  layers. 
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1.  Introduction 

Lithium  ion  cells  are  widely  used  for  portable  equipment 
because  of  their  high  energy  density.  These  cells  employ 
lithium  intercalation  materials  as  their  anode  and  cathode 
active  materials. 

Recently,  there  have  been  several  studies  predicting  the 
cell  voltage  of  cathode  active  materials  in  lithium  ion  cells 
based  on  the  first  principles  quantum  chemistry  [1-4].  The 
internal  energy  was  calculated  using  the  first  principles  quan¬ 
tum  chemistry,  and  the  Gibbs  function  was  obtained.  The 
researchers  in  these  studies  discussed  the  phase  separation 
using  the  Gibbs  function.  We  studied  the  voltage  using  a  more 
simple  method  [5-7],  from  the  Coulomb  potential  created  by 
the  atoms  of  a  cathode  active  material. 

We  used  this  simple  method  to  investigate  the  effect 
of  cathode  particle  size  on  the  cell  potential.  LiNi02  was 
selected  as  a  typical  cathode  of  an  ordered  rock  salt  structure. 

2.  The  voltage  measured  from  Li/Li+ 

We  would  like  to  introduce  our  simple  method  briefly. 
The  following  basic  equations  for  calculating  the  voltage  are 
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derived  from  thermodynamics: 
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where  E  is  the  cell  voltage  (V),  AGr  the  reaction  Gibbs  func¬ 
tion  (J  mol-1 ),  §  the  reaction  coordinate  (mol),  n  the  electrons 
per  molecule  oxidized  or  reduced,  F  the  Faraday  constant 
(C  equiv.-1),  G  the  Gibbs  free  energy  (J),  P  the  pressure  (Pa), 
T the  temperature  (K),  U  the  internal  energy  (J),  V the  volume 
(m3)  and  S  is  the  entropy  (JK-1). 

If  we  consider  Li  insertion  into  the  host  matrix  (HOST) , 
AGr  is  obtained  from  the  difference  of  G  for  the  following 
reaction  divided  by  d$ : 


a£Li  +  Li*(HOST)  ->  Li*^  (HOST) 

From  Eqs.  (l)-(4),  we  have  the  following: 
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The  first  term  in  Eq.  (6)  is  the  internal  energy  term  (Zq nt). 
Let  £/a  be  the  internal  energy  change  caused  by  the  insertion 
of  one  Li.  Then  we  have 
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La  corresponds  to  the  energy  change  that  occurs  during 
the  process  in  which  Li  metal  is  separated  into  one  Li+  atom 
and  one  electron  followed  by  the  insertion  of  one  Li+  atom 
and  one  electron  into  a  host  matrix. 

The  second  term  ( Lw )  is  a  work  performed  against  a  pres¬ 
sure  P.  If  we  consider  the  volume  change  of  the  system,  V  is 
the  volume  of  the  system  and  P  is  the  atmospheric  pressure. 
In  this  case,  the  work  is  very  small.  In  this  paper,  we  assumed 
that  the  matrix  does  not  deform  by  Li  extraction.  Therefore, 
there  is  no  volume  change.  Then,  the  second  term  ( Lw )  can 
be  ignored.  The  lithium  ion  expands  the  host  lattice  when  it 
is  inserted.  Such  a  work  is  included  in  £/a.  We  then  have 
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The  third  term  (Lent)  is  the  entropy  term,  which  is  caused 
by  the  number  of  different  states  of  the  lithium  arrangement 
in  the  host  matrix.  Let  N  be  the  total  number  of  Li  sites 
in  the  host  matrix,  “a”  the  number  of  Li  atoms  inserted  in 
the  host  matrix,  and  NA  Avogadro’s  number.  The  number 
of  ways  that  “a”  indistinguishable  from  Li+  atoms  can  be 
placed  in  N  numbered  sites,  no  more  than  one  to  a  site,  is 
N\/[a !  (N  —  a) !  ] .  If  we  assume  an  equal  probability  of  Li  occu¬ 
pation,  C2p=N\/[a\(N  —  a)\],  we  obtain 


Lent  — 


kT  /ain Q. 


P 


nF  \ 
kTNA 


P,T 


nF 

-In 

RT 


In 


AM 


c N-a -  l)\(a  +  1)! 


AM 


nF 


In 


(N  —  a)\a\ 
N-a 


a  T  1 


RT 

nF 


In 


i-y' 

y 


where  y  =  a/N.  This  is  a  well-known  result  obtained  from 
thermodynamics  [8].  In  this  study,  we  examined  the  entropy 
term  contribution  within  the  framework  of  Eq.  (9).  A  more 
accurate  approximation  of  the  insertion  system  is  discussed 
by  Kudo  and  Hibino  [9] .  Ultimately,  the  voltage  is 
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3.  Phase  separation  and  thermodynamics 

Lrom  Eqs.  (2),  (3),  and  (10),  the  corresponding  molar 
Gibbs  free  energy  [Gm  =  G/(moles  of  matrix)]  is 


(11) 


Fig.  1.  Curve  of  y  lrry  +  (1  —  y)  ln(l  —  _y). 


where  U(y)  is  the  internal  energy  (J  mol-1)  at  a  lithium  occu¬ 
pancy  of  y.  The  curve  of  y  In y  +  (1  —  y)  ln(l  —  y)  was  shown 
in  Lig.  1.  The  tangents  at  y  =  0  and  y=  1  are  +oo  and  —  oo, 
respectively. 

The  relation  between  the  cell  potential,  the  reaction  Gibbs 
function,  and  the  Gibbs  function  are  shown  in  Lig.  2  (left) 
assuming  U(y)  is  a  convex  function.  The  reaction  Gibbs 
function  is  a  mirror  image  of  the  cell  potential.  The  max¬ 
imum  point  A  in  Lig.  2  (left)  (a)  is  the  minimum  point  A 
in  Lig.  2  (left)  (b)  and  the  inflection  point  A  in  Lig.  2  (left) 
(c).  The  minimum  point  B  in  Lig.  2  (left)  (a)  is  the  maxi¬ 
mum  point  B  in  Lig.  2  (left)  (b)  and  the  inflection  point  B 
in  Lig.  2  (left)  (c).  Because  ylny  +  (1  —  y)ln(l  —  y)  versus 
y  curve  shows  a  very  large  tangent  at  y  =  0  and  y  =  1 ,  and 
U(y)  is  a  convex  function,  there  is  a  common  tangential  line 
at  a  lower  side  of  G.  The  points  of  tangency  are  yl  and  y2. 
Let  the  matrix  be  Nm  moles.  When  a  point  y  on  curve  G(y) 
in  Lig.  2  (left)  (c)  is  considered,  it  is  well  known  that  the 
phase  is  separated  to  the  Nm(— y  +  y2)/(y2  —  yl)  moles  of  the 
y  =  y  1  phase  and  the  Nm(y  —  y  1  )/(y2  —  y  1 )  moles  of  the  y  =  y2 
phase.  In  this  case,  the  total  G  of  the  two  phases  is  G  =  Nm 
[(— y  +  y2)G(y  1)  +  (y  -  y  l)G(y)]/(y2  -  y  1).  This  G  value  is  on 
the  contact  line.  The  total  G  of  the  two  phases  is  smaller 
than  the  original  G(y).  Therefore,  the  phase  separation  occurs 
spontaneously.  Lrom  this  consideration,  the  two-phase  region 
is  at  y  1  <  y  <  y2,  and  the  single-phase  regions  are  at  0  <  y  <  y  1 
and  y2<y<l.  Because  the  slope  of  the  common  tangential 
line  is  constant,  the  cell  voltage  is  constant  (Lo)  at  the  two- 
phase  region,  as  shown  in  Lig.  2  (left)  (a).  The  integration  of 
the  reaction  Gibbs  function  gives  the  Gibbs  function.  There¬ 
fore,  the  integrated  value  of  the  reaction  Gibbs  function  from 
yl  to  yl  through  the  original  curve  and  through  the  com¬ 
mon  tangential  line  have  the  same  value.  This  means  that  the 
area  CBL  and  the  area  LAD  are  the  same.  Lrom  this  rela¬ 
tion,  we  can  estimate  the  cell  voltage  (Lo)  at  the  two-phase 
region. 

When  U(y)  is  a  concave  function,  there  is  no  common 
tangential  line  at  the  lower  side  of  G,  as  shown  in  Lig.  2 
(right).  The  insertion  electrode  is  in  a  single  phase  in  this 
case. 


Gm  =  U(y)  +  RT[y  In  y  +  (1  -  y)  ln(l  -  v) ] 
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Fig.  2.  Relation  between  (a)  the  cell  potential,  (b)  the  reaction  Gibbs  function  and  (c)  the  Gibbs  function  [left:  U(y)  is  a  convex  function,  right:  U(y)  is  a  concave 
function.]. 


4.  Phase  separation  of  cathode 


We  studied  the  voltage  using  a  more  simple  method 
[5-7].  We  considered  the  energy  cycle  for  the  reac¬ 
tion  8§Li(crystal)  — >  8§Li+(g)  +  8§e  and  8§Li+(g)  + 
8§e  +  LiANi02  — >  Lix+§^Ni02  [6].  Let  Iu  be  the  ionisa¬ 
tion  potential  (eV)  of  the  Li  metal  (7li  =  —7.107  eV),  and 
@u  be  the  Coulomb  potential  at  a  vacant  Li+  insertion  site. 
The  increase  in  internal  energy  caused  by  the  insertion  of 
an  electron  is  assumed  to  be  —  +  In\  ,  where  <7>Ni  is  the 

Coulomb  potential  of  the  electron  insertion  site  and  /ni 
is  the  LUMO  energy  of  an  isolated  Ni  atom  in  a  vacuum 
(7]sji  —  — 29. 106  eV).  When  <7>Ni  is  calculated  for  a  cathode 
before  the  Li+  insertion,  we  have  to  add  the  contribution  of 
the  inserted  Li+.  Thus,  we  have  to  add  —  14.399/TVi-Ni  to 
0Ni>  where  Ru-Ni  is  the  distance  between  the  inserted  Li+ 
and  the  electron.  Finally,  we  have 
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where  the  sign  after  7pi  is  +  if  there  is  already  Li  at  the  center 
of  k  (extraction  of  Li)  and  —  if  there  is  no  Li  at  the  center  of 
k  (insertion  of  Li).  With  Li-ion  cells,  n=  1. 

Let  rij(A)  be  the  distance  from  a  certain  atom  k  to  the 
surrounding  atoms  of  species  i  in  an  ionic  crystal,  and  Z*  be 


the  charge  number  of  species  i.  Then  the  Coulomb  potential 
at  atom  k  is  expressed  by  the  following  equation: 


0k  =  14.399 


(13) 


On  the  surface  of  the  LiNi02,  “Z;” s  at  plane,  edge,  and 
corner  are  assumed  to  be  1/2Z;,  1/4 Z;  and  1/8Z;,  respectively. 
<7>k  was  calculated  using  VIP  BASIC  (Mainstay)  for  Power 
Macintosh.  The  value  of  r#  was  obtained  from  the  crystal 
structure  of  a  cathode  active  material.  The  summation  of  r# 
was  undertaken  for  all  the  ions  until  the  value  converged. 


5.  Voltage  prediction  of  LiNi02 

We  calculated  the  voltage  of  LiNi02  [6].  The  difference 
between  the  ideal  structure  used  for  our  calculations  and  the 
real  structure  is  shown  in  Table  1 .  We  assumed  in-plane  order¬ 
ings  of  lithium  ions  as  shown  in  Fig.  3.  In  each  area,  Li+-ions 
are  assumed  to  be  distributed  at  random.  This  means  that 
all  the  sites  are  occupied  by  Liv+  and  ev_ ,  where  y  is  the  site 
occupancy  in  each  area.  Then,  for  7qnt,  we  have  a  straight  line 
connecting  two  voltages  (7qnt)  at  each  end  point  of  the  area. 
The  Gibbs  function  and  the  cell  voltage  of  LiNi02  are  shown 
in  Figs.  4  and  5.  The  experimental  results  of  LiANi02  are  two 
phases  (R-3m)  at  0<v<0.25  with  4.2  V,  a  single  phase  (R- 
3 ni)  at  0.25  <v< 0.45  with  4. 2-3. 8  V,  a  single  phase  (C2 hn) 
at  0.45  < v< 0.75  with  3. 8-3. 7  V,  and  a  single  phase  (R-3m/) 
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Table  1 

Differences  between  the  ideal  structure  and  the  real  structure 

Nearest  length  (A) 


Ni-O 

Li-0 

Ni-Ni 

Li-Li 

O  layer  distance 

LiNi02 

Ideal  structure 

\/2  r 

\/2  r 

2r 

2  r 

(4/s/6 )r 

():r  =  V 2  A 

(2.000) 

(2.000) 

(2.828) 

(2.828) 

(2.309) 

Real  structure 

1.957 

2.126 

2.875 

2.875 

2.36 

0  x  0.25 


0.25  x  0.50 


0.50  x  0.75 


0.75  x  1.00 


O  :  Vacant  Li  site  %  •  Occupied  Li  site 


Fig.  3.  In-plane  orderings  of  lithium  ions. 


at  0.75  <x<  1  with  3. 7-3. 5  V.  Our  calculated  phase  separa¬ 
tion  was  in  agreement  with  the  experiments  at  0<x<0.25. 
This  bad  result  probably  occurred  because  we  assumed  that 
the  matrix  does  not  deform  by  the  insertion  of  lithium  ions 
and  electrons  (iua  =  0).  The  calculated  voltage  at  1  was 
in  very  good  agreement  with  the  experimental  value,  because 
the  deformation  of  the  matrix  was  very  small  at  the  beginning 
of  the  extraction  of  Li  from  the  LiNi02  matrix. 


6.  Voltage  prediction  of  nano-size  LiNi02 

0k  was  calculated  using  Eq.  (13).  Usually,  the  summation 
of  rij  was  undertaken  for  all  the  ions  until  the  value  converged, 
but  when  LiNi02  is  very  small,  there  are  not  enough  ions 
around  the  atom  k.  Therefore,  the  0k  becomes  a  different 
value  of  the  converged  value.  The  0k  depends  on  the  position 


in  LiNi02  particle.  In  this  paper,  0k  is  calculated  at  the  center 
of  LiNi02  particle.  The  calculation  of  voltage  was  performed 
for  E[nt  of  LiNi02  and  Lio.75Ni02.  The  two-phase  voltage 
fromx=  1  to  0.75  (FiANiC>2)  was  also  calculated.  Those  three 
voltages  are  shown  in  Fig.  5.  As  already  mentioned,  the  most 
reliable  voltage  is  Emi  of  LiNi02.  To  increase  the  size  of 
LiNi02  from  a  unit  cell,  we  examined  three  different  cases. 

6.1.  Increasing  the  size  ofLiNi02  while  maintaining  a 
similar  shape  of  the  unit  cell 

Shape  of  LiNi02  cathode  was  shown  in  Fig.  6,  when 
a  similar  shape  of  the  unit  cell  was  maintained.  The  volt¬ 
age  increased  as  the  size  of  FiNi02  decreased,  as  shown  in 
Fig.  7.  In  this  calculation,  the  outside  layers  are  Fi  layers. 
Emi  of  FiNi02  increased  very  rapidly  from  4.3  to  12.3  V  as 
the  size  of  FiNiC>2  decreased.  Emt  of  FiNiC>2  was  4.3254  V 


Fig.  4.  Gibbs  function  of  LixNi02. 
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Fig.  5.  Cell  voltage  of  LiANiC>2.  Experimental  data  is  from  Ohzuku  [10]. 

at  b  =  139  A  and  4.3366  V  at  b  =  97  A.  This  means  that  Emi 
of  LiNiC>2  increased  for  11.2  mV  by  the  decrease  of  the  size 
from  b  =  139  to  97  A.  It  is  not  very  difficult  to  measure  the 
10  mV  difference.  Therefore,  if  we  synthesize  100  A  (10  fxm) 
particles  of  LiNi02  experimentally,  we  can  check  the  volt¬ 
age  change.  We  are  now  conducting  experiments  to  prove  our 
theory. 

The  case  of  the  outside  oxygen  layers  was  also  calculated 
(Figs.  8  and  9),  because  Li  must  be  inserted  between  two 
oxygen  layers.  The  result  was  similar  to  the  case  of  the  Li 
outside  layers.  Emi  of  LiNi02  was  3.6598Vat/?=139A  and 
3.6739  V  at  b  =  97  A.  The  difference  in  voltage  is  14.1  mV. 


Fig.  6.  Shape  of  LiNi02  cathode  when  a  similar  shape  of  the  unit  cell  was 
maintained. 


Fig.  7.  Zsjnt  of  FiNi02  (Li  layer  is  outside  maintaining  a  similar  shape  of  the 
unit  cell). 


0  50  100  150 


Length  of  b  (A) 

Fig.  8.  Ejnt  of  FiNi02  (oxygen  layer  is  outside  maintaining  a  similar  shape 
of  the  unit  cell). 


Length  of  b  (A) 

Fig.  9.  Fjnt  at  Fi  =  1  of  FiNi02  (oxygen  layer  is  outside  maintaining  a  similar 
shape  of  the  unit  cell). 


6.2.  Increasing  the  size  ofLiNiC>2  while  changing  the 
in-plane  area  and  keeping  the  number  of  layers  (20 
layers) 

The  voltage  increased  as  the  size  of  LiNiC>2  decreased,  as 
shown  in  Fig.  10.  In  this  calculation,  the  outside  layers  are 
Li  layers.  Emi  of  LiNi02  increased  very  rapidly  from  4.6  to 
12.6V  as  the  size  of  LiNiC>2  decreased. 
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0  10  20  30  40  50  60 


Length  of  a  (A) 

Fig.  10.  E-mt  of  LiNi02  (Li  layer  is  outside  changing  the  in-plane  area  and 
keeping  the  number  of  layers). 


Length  of  b  (A) 


Fig.  11.  E-mi  of  LiNi02  (Li  layer  is  outside  maintaining  the  in-plane  area). 


Length  of  b  (A) 

Fig.  12.  E-m{  of  LiNi02  (oxygen  layer  is  outside,  maintaining  the  in-plane 
area). 


6.3.  Increasing  the  size  of  LiNiO 2  while  maintaining  the 
in-plane  area  (20  x  20  unit  cells )  and  changing  the 
number  of  layers 

The  voltage  decreased  slightly  as  the  size  of  LiNiC>2 
decreased,  as  shown  in  Fig.  11.  In  this  calculation,  the  out¬ 
side  layers  are  Li  layers.  Emi  of  LiNi02  decreased  from  4.3 
to  4. 1  V  as  the  size  of  LiNiC>2  decreased. 

We  also  calculated  the  case  of  outside  oxygen  layers 
(Fig.  12).  E-mi  of  LiNiC>2  decreased  very  rapidly  from  4.3 
to  —12.6V  as  the  size  of  LiNiC>2  decreased. 

From  a  practical  point  of  view,  the  decrease  of  the  voltage 
is  not  good.  We  conclude,  therefore,  that  reducing  the  number 
of  layers  is  not  a  good  way  to  synthesize  nano-sized  LiNi02. 
However,  we  have  to  find  a  new  electrolyte  which  have  a 
higher  oxidation  potential,  when  the  high  voltage  cathode  is 
used. 

7.  Conclusion 

Based  on  the  results  of  this  study,  we  conclude  that  the 
voltage  changes  when  the  size  of  LiNiC>2  becomes  less  than 
about  10  nm.  The  voltage  increased  as  we  decreased  the  in¬ 
plane  size  and  decreased  as  we  decreased  the  number  of 
layers.  However,  depends  on  the  position  in  LiNiC>2  par¬ 
ticle.  And  the  matrix  may  easily  deform  during  Li  extraction 
when  the  size  of  LiNiC>2  particle  becomes  very  small.  We  are 
now  studying  the  contribution  of  these  factors. 
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